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Abstract: The fragmentation of a light parton i to a jet containing a hght energetic hadron 
/i, where the momentum fraction of this hadron as weh as the invariant mass of the jet is 
measured, is described by "fragmenting jet functions" . We calculate the one-loop matching 
coefficients Jij that relate the fragmenting jet functions to the standard, unpolarized 
fragmentation functions for quark and gluon jets. We perform this calculation using 
various IR regulators and show explicitly how the IR divergences cancel in the matching. We 
derive the relationship between the coefficients Jij and the quark and gluon jet functions. This 
provides a cross-check of our results. As an application we study the process e'^e~ — )• X-n'^ 
on the T(4S') resonance where we measure the momentum fraction of the vr"'' and restrict 
to the dijet limit by imposing a cut on thrust T. In our analysis we sum the logarithms of 
r = 1 — r in the cross section to next-to-next-to- leading- logarithmic accuracy (NNLL). We 
find that including contributions up to NNLL (or NLO) can have a large impact on extracting 
fragmentation functions from e^e~ — )• dijet + h. 
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1 Introduction 

In single-inclusive hadron production, an energetic parton i = {g,u,u,d, . . .} produces an 
observed energetic hadron h and accompanying hadrons X. Factorization theorems allow one 
to identify perturbative (calculable) and non-perturbative (universal) contributions to these 
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processes. For example, in e~^e^ — t- Xh at a high center-of-mass (cm.) energy Q it has been 
proven that, to ah orders in and at leading power of Aqqd/Q, the cross-section has the 
following factorized form (see e.g. ref. [1]) 

g = ao Y: l''^a{Q/-,,)DHx,,). (1.1) 

i=g,u,u,d,... 

Here, z = 2Eh/Q is the energy fraction of the hadron h in the cm. frame, co is the Born 
cross section, and // is the MS renormalization scale. The coefficient functions Q incorporate 
the short-distance partonic process producing the fragmenting parton i: they are calculable 
in perturbation theory and independent of the observed hadron h. The long-distance physics 
of the hadronization resides in the non-perturbative fragmentation function D^{x,^), which 
is the number density of hadrons of type h in the "decay" products of the parton i, for a 
specific value of x [2-5]. The convolution variable x in eq. (1.1) is the fraction of the energy 
of the parent parton i carried by the observed hadron h. At leading order (LO) the hard 
partonic process is e'^e~ — )• qq and x = z. Beyond LO, radiation will be emitted before the 
parton i produces h, and thus x > z. 

In experimental studies of fragmentation, additional measurements on the hadronic final 
state X may be needed. For example, the Belle collaboration studies light-quark fragmenta- 
tion by restricting to dijet final-state configurations, which removes i?-meson events from the 
data sample on the T{AS) resonance [6]. This is achieved by imposing a cut on thrust, which 
is an event shape variable defined as [7] 

T = max,~^iM, (1.2) 

\Pi\ 

where the sum is over all final-state particles. In terms of r = 1 — T, a more convenient 
quantity to describe dijet events, r close to corresponds to configurations with two narrow, 
pencil-like, back-to-back jets; while the other extreme r = 1/2 corresponds to a spherically 
symmetric event. Since at ^/s = 10.58 GeV the B mesons decay nearly at rest in the cm., a 
thrust cut of T < 0.2 removes 98% of the B data leaving the thrust distribution dominated 
by the fragmentation of light (uds) and charmed quark pairs [6]. This cut on r constrains 
the (squared) invariant masses Si of the final-state jets. Indeed, in the dijet limit, 

Sa + Sh k 

where k is the contribution from soft radiation between jets. 

We focus on such restrictions on the hadronic final state, by studying the (spin-averaged) 
fragmentation of a light hadron h inside a collimated jet originating from a light parton i, 
when the jet invariant mass is constrained. These features cannot be solely described by 
D^[x,fi), which only depends on the momentum fraction x: in our case fragmentation is 
probed at a more differential level. In ref. [8] a novel "fragmenting jet function" Q^{s,z,ii) 
was introduced which depends both on the fragmentation variable z and on the invariant 
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mass s of the collinear radiation that forms the jet. The relevant hierarchy of scales is given 
by m/j <^ t/s <^ Ejet, where E'jet is the jet energy. This hierarchy allows us to employ Soft- 
Collinear Effective Theory (SCET) [9-12], which is an effective field theory of QCD suitable for 
processes with well-separated energetic jets. Collinear and soft degrees of freedom describe, 
respectively, the energetic radiation inside jets and the soft emissions between them. A 
different collinear sector is associated with each jet. The collinear sectors and the soft sector 
decouple at leading power [12]. This leads to factorization formulae for inclusive observables 
at high energies which are characterized by convolutions of jet functions J{si,fi) (describing 
the invariant mass distribution of each jet) with a soft function S encoding the contribution 
of the soft degrees of freedom, see e.g. eq. (1.6) below. 

The fragmenting jet function (s, z, n) has features of the standard fragmentation func- 
tion D^{x, //) and the leading inclusive jet function Ji(s, ^), which is calculable in perturbation 
theory. As was shown in ref. [8] , the following simple replacement rule holds 

•Ji{s. /") ^ GHs. fi) dz , (1.4) 

which allows us to obtain factorization formulae for semi-inclusive processes with fragmenta- 
tion within a jet, from the corresponding inclusive ones. The factor 2(27r)^ is related to the 
normalization of and to the phase space factor for the hadron h. 

Here we will focus on the relation between Q^{s,z,fj,) and D^{z,^). At leading order 
in Aq^j-,/s <^ 1, the fragmenting jet function can be expressed as a convolution between 
short-distance coefficients and the standard fragmentation functions at the scale ~ y/s [8], 




The Jij describe the emission of collinear radiation, forming a jet with invariant mass s, 
within which the non-perturbative, long-distance fragmentation process takes place. 

In this paper we present the one- loop calculation of the matching coefficients J'ij , where 
the initiating parton i can be either an (anti)quark or a gluon. This completes the picture 
detailed in ref. [8] with the information necessary to relate factorization theorems for semi- 
inclusive processes, where the jet invariant mass is probed, with the standard -D^(x,/i). 

In the presentation of our results particular attention will be devoted to show how the 
infrared (IR) divergences of the partonic Gi and Di cancel in the matching. Cross-checks of our 
results are provided by the known anomalous dimensions of the fragmentation functions and 
the (fragmenting) jet functions, as well as by a relationship between Qi(s,z,fj,) and Ji{s,fi), 
which we work out in sec. 2.5. As we will see explicitly, the J7ij contain double logarithms, 
i.e. contributions of the form a" ln'"(s//ij) with m < 2n, so that eq. (1.5) should be evaluated 
at ~ -y/s to avoid the breakdown of the standard perturbative expansion. 

We will illustrate our results through a numerical analysis of the process e~^e~ — )• Xh 
where we restrict to the dijet limit by a cut on r, as in the Belle study of light quark 
fragmentation mentioned above. Using eq. (1.4) in the leading-order factorization theorem 
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Figure 1. Schematic display of the factorization in eq. (1.6) for the fragmentation process e+e" — >■ 
dijet + /i. The cross denotes the short-distance process e+e~ — >• qq producing the back-to-back jets. In 
one of these jets a hadron of type h is observed and its momentum fraction z is measured. Each event 
may contribute more than once to the cross-section if the final state contains several of these hadrons 
h{zi), . . . , h{zn)- The dijet limit restricts the radiation to be either collinear or soft, drawn respectively 
in black and red color. At leading power, the collinear radiation is described by a (fragmenting) jet 
function, and the two jets only interact through soft radiation, described by the soft function. 



for the cross-section where an inclusive measurement of thrust is performed for r ^ 1 [13-16], 
we obtain: 



drdz ^2(27r)3 



H{Q ,/x) / dsadsbdk 



Sg + Sb 



k 
Q 



On 



2(27r): 



l + 0(r) 
dx 



Jc 



qj 



X D'l{x,fi)QSr[QT 



Sa + Sb 



Q 



A2 



^) Jq{Sb,fJ-) + Jq{Sa,^j) Jqj\Sb 



Z 
X 



(1.6) 



see fig. 1. Here eq. (1.3) is incorporated through the (^-function. The first line receives power 
corrections of 0{t) and in the second line we also have 0[Aqqy)/{tQ'^)] corrections from using 
eq. (1.5). We will only consider the contribution from the light quark flavors q = u,u, d, d, s, s. 
The gluon fragmenting jet function does not appear in eq. (1.6), but the gluon fragmentation 
function does contribute because the sum over j includes j = g. The normalization factor 
(Tq is the tree- level cross-section for the electroweak process e^e~ — )• (7 , Z) — )• given in 
eq. (C.l), which depends on the quark flavor. Since we assume that it is not known whether 
the observed hadron h fragmented from the quark or the antiquark initiated jet, we have a 
sum over both possibilities in the factorization theorem. 

In eq. (1.6), the hard function H{Q'^,ii) encodes virtual effects arising from the production 
of the qq pair at the hard scale — Q, and is given by the square of Wilson coefficients 
in the matching of the relevant QCD onto SCET currents. The (real and virtual) collinear 
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radiation of the jet from which the hadron fragments is described by or Jij whereas 
the jet in the opposite hemisphere is represented through an inclusive jet function. The 
characteristic scale of these functions is the jet scale /ij ~ \/tQ. Finally, the soft function 
Sr{k, fj.) describes the contribution to the hemisphere masses (and therefore to thrust) due 
to soft parton emissions. Sr is defined through the vacuum matrix element of eikonal Wilson 
lines and the corresponding soft scale is ^5 ~ tQ. 

In the two-jet limit r <C 1, the cross section in eq. (1.6) contains large double logarithms 
a" In™ r (m < 2n) , which need to be resummed to make reliable predictions and uncertainty 
estimates. In our effective field theory approach this is achieved by evaluating the hard, 
(fragmenting) jet and soft functions at their natural scales ^h, and ns respectively, where 
they contain no large logarithms, and by running them to a common scale fi using their 
respective renormalization group equations (RGEs). In eq. (1.6) all functions except the 
matching coefficients Jij have already been studied in the literature. Our calculation therefore 
provides the missing ingredient necessary to sum these logarithms up to next-to-next-to- 
leading-logarithmic (NNLL) accuracy. We will discuss resummation effects in eq. (1.6) for 
the case of unpolarized single charged pion production when the soft scale ns is perturbative. 

We stress that our results have a broad range of applicability. For example, eq. (1.4) 
can be directly utilized when the hadronic final state is characterized via the event shape A^- 
jettiness tn [17], since in this case the jets are described by the standard jet functions. This 
variable measures how A^-jet-like an event is, and can be used to veto unwanted additional 
jets by requiring r^v < t™* <^ 1, which is the region of validity of the A^-jettiness factorization 
theorem. In a more exclusive approach where jet algorithms are used, jet-algorithm-dependent 
jet functions arise [18, 19]. However, even for A^-jettiness one may use a jet algorithm to 
determine the jet energies and directions, since all reasonable jet algorithms agree on these 
quantities in the exclusive A^-jet regime tat <C 1, up to power corrections of 0{Tj\f) [17]. 

The paper is organized as follows. In section 2 we set up the theoretical framework, review 
the SCET definitions of fragmenting jet functions and standard fragmentation functions, and 
discuss their relationship and renormalization properties. We also discuss the relationship 
between the fragmenting jet function and jet function in detail. Our results for the matching 
coefficients Jij are given in section 2.4. In section 3 we present our calculation of the quark 
matching coefficients Jqi at one loop, where we show how the IR divergences cancel in the 
matching procedure. Section 4 is devoted to the gluon case at next-to-leading order (NLO). 
Numerical results for the fragmenting jet functions using eq. (1.5), as well as, a numerical 
analysis of pion fragmentation in dijet-like e~^e~ — )• Xtt~^ with a cut on thrust are contained 
in section 5. Conclusions and outlook are given in section 6. Useful mathematical identities 
are given in appendix A. An alternative calculation of J7gj using the optical theorem and 
a different IR regulator can be found in appendix B. All the ingredients necessary for our 
numerical analysis are collected in appendix C. 
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2 Fragmentation within an Identified Jet 



In this section we start by setting up the theoretical framework of our analysis, and introduce 
the SCET ingredients relevant for this paper. We give the definitions of quark and gluon 
fragmentation functions in SCET and discuss their renormalization. We then focus on quark 
and gluon fragmenting jet functions, discuss their renormalization and their relationship with 
the standard fragmentation functions. At the end of this section, we consider their relationship 
with the jet function, which provides a powerful cross-check on our calculations. 

2.1 SCET Ingredients 

Light quark or gluon fragmentation within an identified jet is governed by three different 
scales: the (perturbative) hard scale set by the jet energy Exh, the intermediate (perturba- 
tive) jet scale given by the jet invariant mass mxh, and the soft scale of order ~ m\y^/ Exh, 
with the hierarchy mh <C mxh Exh] here we will always consider the hadron mass rrih to 
be negligible as in ref. [8]. After integrating out the hard dynamics that initiates the jet, we 
are left with collinear and soft modes in the Xh system. Therefore, SCET [9-12] - which is 
an effective field theory of QCD that describes physics of collinear and soft degrees of freedom 
- is well suited for this analysis. 

Since the invariant mass of a jet is much smaller than its energy, the jet constituents are 
collimated and conveniently described using light-cone coordinates. To this end we introduce 
a light-cone vector n'^ whose spatial part is along the jet axis, and another light-cone vector 
such that v? = = {) and n-n = 2. Any four- vector can then be decomposed as 
= {p'^,p~ ,Pj_) with p^ = n-p, p~ = n-p and pj_, which contains the components of p^ 
perpendicular to n^^ and n^. The momentum p^^ of a particle within the jet scales collinear ly, 
i.e. p^ = ~ p^(A^,l,A), where A ~ rnxh/Exh ^ 1 is the SCET expansion 

parameter. For the soft degrees of freedom, the momentum scales like = {q^,q^,qj_) ~ 



The collinear momentum is separated into a large part and a small residual part 



with = i0,pj,pe±) ~ pj{0, 1, A) and p^ = ipt,Pr,Pr±) ~ ^^)- The SCET fields 



for n-collinear quarks and gluons, Cn,piiy) and An^p^di) respectively, are labeled by n and 
the label momentum p£. Their argument y is conjugate to the small residual momenta. A 
derivative acting on these fields picks out the residual momentum dependence, id^ ~ ~ 
X'^pJ , while label momentum operators Vn = n ■ Vn CPnx ) return the sum of the minus 
(perpendicular) label components of all n-collinear fields on which they act. 

Interactions between collinear fields cannot change the direction n but change the mo- 
mentum labels. It is therefore convenient to use the short-hand notation 



P-(A2,A2,A2). 



n-piY +Pe±+Pr 



(2.1) 





(2.2) 
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In the sum we explicitly exclude the case = to avoid double-counting of the soft degrees 
of freedom (which are described by separate soft quark and gluon fields). In practice, when 
calculating matrix elements, this is implemented using zero-bin subtractions [20] or alterna- 
tively by dividing out matrix elements of Wilson lines [21-23]. We will study the zero-bin 
subtractions in detail since they play an important role in our calculation. This is most ex- 
plicitly seen in section 3, where we use a gluon mass and a (^-regulator [24] to regulate the IR 
divergences in the one-loop quark fragmenting jet function. 

Collinear operators are built out of products of fields and Wilson lines that are invariant 
under collinear gauge transformations [10, 11]. The basic building blocks are the collinearly 
gauge-invariant quark and gluon fields, defined as 



Xn[y) = wl{y)Uy) 



B'^^{y) = -^[wl{y)\D^^^Wn{y)], 



(2.3) 



where i-D^_|_ = + g^^i _L-collinear covariant derivative. The collinear Wilson line 



(2.4) 



Wniy) = Yl ^^p(" J" ""-My) , 

'-perms " 

sums up arbitrary emissions of n-collinear gluons from an n-collinear quark or gluon, which 
are 0(1) in the power counting. 

At leading order in the SCET power expansion, the interactions of soft gluons with 
collinear fields exponentiate to form eikonal Wilson lines. The soft gluons can thus be decou- 
pled via the BPS field redefinition [12] 



xffL(y) = ^rI(y)xn,c.(y), 
K''^i.iy) = yli.y)Kui.iy)yn{y). 



(2.5) 



The collinear fields we consider in this paper are those after this decoupling, and we drop the 
superscript (0) for notational convenience. Here Yn{y) is a soft Wilson line in the fundamental 
representation 



Yn{y) = Pexp 



\g I dun-Ausiy + un) 





9 



perms 



ire 



d 



n-Aus{y) 



(2.6) 



The symbol P in eq. (2.6) denotes anti-path ordering of the color generators along the in- 
tegration path. On the second line we write the Wilson line in momentum space akin to 
eq. (2.4). 

2.2 Pragmentation Functions 

The fragmentation functions D'^{x) characterize the factorization theorems that describe 
high-energy single-inclusive hadron production processes at leading power [1] , where no prop- 
erties of the jet are probed, see e.g. eq. (1.1). These functions encode the non-perturbative 
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information on how the energetic parton i (either a gluon or an (anti)quark of a certain fla- 
vor) produces the observed hadron h which carries a fraction x of the initial parton's large 
light-cone momentum component. 

Let and denote the parton and hadron momenta, respectively. In a frame where 
k± = 0, the hadron has pf^ = xk~ and = {p^x + iTT-h) /Ph ■ With the gauge choice n- A = 0, 
the bare unpolarized quark fragmentation function has the following operator definition in 
QCD [5] 

D'^Mrei^) = If d'pi |^|^^e*^-^^^/2^^tr[|(0|^(y+,0,y^)|XM(X/.|V;(0)|0) 

(2.7) 

where ip is the quark field quantized on y~ = and the trace is taken over color and Dirac 
indices. The factor l/{2Nc), where A'c = 3 is the number of colors, comes from averaging 
over the color and spin of the parent parton. The state \Xh) = \Xh{pfi)) contains a hadron 
h with momentum ph, and a sum over the polarizations of h is assumed. Boost invariance 
along the non-_L direction implies that D can only be a function of x = pj^ /k~ and not pj^ or 
k~ individually. According to factorization at leading power, the sum over the accompanying 
hadrons X is dominated by jet-like configurations for the \Xh) states [1]. 

In SCET notation, the fragmentation function takes on the following form [8] 

^',bare(^) = ^ / E [f -^fe v)'^' (P^r) (0| '^CPx Xn(0)] | X/l> Xn (0) 1 0> 

(2.8) 



Here, Xn is the n-collinear quark field in eq. (2.3) that contains a Wilson line, making this 
definition (collinearly) gauge invariant. The V and V± operators pick out the 0{X^) and 
0{X) label momentum of the field, while the continuous 0{X'^) residual components of the 
jet momentum are denoted by Pxhr- notation = + Ph- 

The QCD definition for the bare gluon fragmentation function, in the n ■ A = gauge 
and in a frame where pj^f^ = 0, for d space-time dimensions, reads [5] 

1 f.2i f dy+ d'^y± ik-y+/2 



^sMrei^) (d_2)(iV2_i)p-y ^^^y 2(27r)3 

xY,n^n''{0\G^,^{y+,0,y^)\Xh){Xh\Gt'''{Om (2.9) 



X 



where G^i, = C^^, T"" is the QCD field-strength tensor and an average over colors and the 
{d — 2) polarizations of the gluon is performed. The corresponding expression in SCET is 
given in gauge invariant form by 

^3,bare(^) = ' - 2){N^, - I) X j "^^^"^ ^ '^^^^'^'^^ ^^P^Kr) 

X <0|[5,,^5o,P,^::f(0)]|X/i)(X/i|^^_^(0)|0> . (2.10) 
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The operator products in the definitions of D^{x) and Dg{x) are singular and require 



renormahzation. The renormahzed fragmentation functions are defined through [5] 

dx' 

j=g,u,u,d,. 



Dt^.re{x)= E r^^^(5'/^)^i(^''^)- (2-11) 



Throughout this paper [i denotes the scale of dimensional regularization in the MS scheme. 
The renormalization-group equation (RGE) for D^{x,^) follows from eq. (2.11) 

%A^(-'^) = Ef (2-12) 
with anomalous dimension 

Here, the inverse of the renormahzation factor {Z^)^f} is defined through 

k -'^ 

In our perturbative calculations we will replace the hadron h by either a quark or a gluon. 
The other ingredients of a factorization theorem are not affected by this, e.g. the hard, jet and 
soft function in eq. (1.6) are the same in both cases. Since the cross section is an observable, 
the /U dependence must cancel between all the factors of a factorization theorem. From this it 
follows that the renormahzation and anomalous dimension of D are the same if /i is a hadron 
or a parton. 

We will use the variable p for the momentum of the parton that replaces the hadron h. 
Denoting the discrete label parts of the momentum p^^ by p'^ and pii±_ and the continuous 
residual parts by pr, the partonic fragmentation functions are at tree- level given by 

Df^i^) = IT.I dV^tr[|5^^^-5o,^x5(p-)52(^±) (0|en(0)|g„(p)>(g„(p)|en(0)|0> 

1 _ Vf/i 



2x 



5(^^-p-)tr[|^<(pX(p)] =5{l-x) 



X {0\A';^'l{0)\9n{p)){gn{p)\A^n±,M0) 



6{u-p-) Y.elAp)-^n±ip) = 5(1 - x) . (2.15) 



(d-2)x 

^ ' pol 

Here we recombined residuals and labels into the continuous p^ , via 



^.,p-^iPr)=Si^-p-) ■ (2.16) 



- 9 - 



We will use this relation in the rest of the paper. The partonic Dq and Dg vanish at tree 
level. At one-loop their UV divergences lead to the mixing of quark- and gluon fragmentation 
functions. The details of this calculation are given in section 3. 

The one- loop renormalization factors zl^^^\x,fi) can be easily extracted once the one- 
loop partonic result is known. Expanding the partonic version of eq. (2.11) to one-loop, 



E 

k 



X' 

dx^ 



X 



D(l) 



„^,)Df\x',^,) 



Sik S I 1 



X' 



Dpix',,) + Z-^^^ 



\ 6kj5{l - x') 



Dl^'\x,^) + zf/'\x,t,), 



(2.17) 



where the superscripts (0) and (1) denote the tree- level and one- loop expressions, respectively. 
From eqs. (2.13) and (2.17) the one-loop 7j^(x,/x) can then be obtained straightforwardly. 
At 0{as) the space- like and time-like Altarelli-Parisi evolution kernels are related to each 
other via a simple analytic continuation rule [25, 26] based on symmetries of the relevant 
diagrams under crossing, and via the so-called Gribov-Lipatov reciprocity relation which 
connects space-like and time-like structure functions in their respective physical regions [27, 
28]. As a consequence, the 7j^(x,/x) coincide with the anomalous dimensions of the parton 
distribution functions at one loop (where the roles of incoming and outgoing partons are 
interchanged). Therefore, 



^qg{X, 



TT 



e{x)Pggix) , 

e{x)Pgq{x), 



1 



e{x)\CAPgg{x) + -^0 6il 

as{n)TF 



TT 



e{x)Pqg{x) 



(2.18) 



where /3o = (HCa — 4n/Tp)/3, is the lowest order coefficient of the QCD /3-function. The 
splitting functions are [29] 



1 



Pgq{x) = e{l - X) 



= (l + x2)£o(l 
l + {l-xf 



1 



+ x{l -x) 



Pgg{x) = 2x£o(l -x) + 2e{l- x) 

Pqg{x) = e{l-x)[x^ + {l-xf], (2.19) 

in which we do not include the usual color factors for future convenience. The plus-distribution 

fl-x) 



1 



(2.20) 
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is defined in eq. (A. 2). In section 3 we will show how our partonic calculation of D^q^\x,^) 
also leads to eq. (2.18). The relation 7j^(x, ^) = "fj^ix, fj.) was shown no longer to hold beyond 
one loop in dimensional regularization and the MS scheme in ref. [4], see also the discussion 
in ref. [30]. 

2.3 Fragmenting Jet Functions 

A high-energy single-inclusive hadron production process, where the invariant mass of the 
jet initiated by a light parton is measured, is described by a fragmenting jet function. We 
consider the case that the hadron mass is negligible (compared to the jet mass) and work in a 
frame where the perpendicular momentum of the jet is vanishing. For a fragmenting quark [8] 



x{0\K,v^0Mn{y-)]\Xh){Xh\xnm\0)\ . (2.21) 

In the second equality we performed a translation of the collinear field Xniu), whose argu- 
ment is associated with residual momenta, and carried out the integrals in and y*". The 
integration over y~ fixes the partonic jet invariant mass s = k~^uj. Analogously, in the case 
of a gluon-initiated jet, in d space-time dimensions, 

X (0|[i„,pVj'xS;i°(!/-)l|X/.>(Xft|S;j_./0)|0>. (2.22) 

Evaluating the partonic fragmenting jet functions at tree- level, we find 



gf\s, z) = 5{k-^) 5{u - p-) 1 tr [| ^ un{p)un{p)\ = 2{2^f5{s)5{l - z 



spins 

^3, 



gf)^s, z) = s{u - p~) ^ = 2{27:f5{s)5{l - z) . (2.23) 

^ pol 

Here we have used the fact that p"*" = 0, due to the on-shell condition = and the choice 
of frame, which sets = 0. 

A consequence of eq. (1.4) is that the renormalization and RG evolution of these two 
functions are the same 

GiMre{s,z)= rds'Zlj{s-s',f,)g^{s',z,fi), (s, /.) = Z}(s, /x) , (2.24) 

Jo 

where the index i is not summed over. In particular, the renormalization of does not affect 
its z-dependence and does not mix quark and gluon fragmenting jet functions, at any order 
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in perturbation theory. We will see this explicitly in our one-loop calculation in sections 3 
and 4. The corresponding RGE is given by 



where 



and (Z'r 



' ■ 1 d • 

ds' {Zg) {s - s' , n) fi—Zg{s' , fl) , 



is defined as 



/V {Zl,)-\s-s',^^) Ziyis', f,) = 5{s). 
Jo 



The structure of the anomalous dimension of the jet function jj{as) implies 



+ fg{as)S{s), fg{as) = fj{as) , 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



where the plus distribution Co is defined in eq. (A.2). The cusp anomalous dimension 
r*^gp(as) [31] and the non-cusp part of the anomalous dimension 7g(as) are collected in 
the app. C, to make the paper self-contained in view of the numerical analysis in sec. 5. A 
cross-check of our partonic one-loop calculation of Qj^^\s, z, fj,) will be provided by the one- 
loop evolution kernels 7g(s,/u) in eq. (2.28). From the anomalous dimension we can obtain 
the one-loop renormalization factor Zg^^ through eq. (2.26). This can be compared with our 
calculation by using eq. (2.25) expanded to one loop: 



ds' 
ds' 



) (. - /x) gj''> is', z, /.) + zl''> {s - s', ^^) gr is', z, ^^) 



6{s - s') gf^\s', z, fi) + )(s - s', fi) 6ij 6{s') 6{l - z) 



gi''>{s,z,,i) + Z'g''>{s,ij)5.,j5{l-z). 



(2.29) 



Here we have noted that gl contributes at tree level only when i = j. Therefore for i ^ j, gj 
are UV finite at one loop. 

2.4 Results for Matching onto Fragmentation Functions 

By performing operator product expansions of the fragmenting jet functions in eqs. (2.21) 
and (2.22) about the y~ — )• limit, we can match onto the low-energy matrix elements in 
eqs. (2.8) and (2.10) that correspond to the fragmentation functions. This amounts to the 
SCETi onto SCETn matching (illustrated in fig. 2) at the intermediate scale provided by the 
jet invariant mass fij ~ y/s: 



j=g,u,u,d,. 



1 dz' 



Ji 



'-,f,j)D^;{z',f,j) 



l + O 



^^QCD 



(2.30) 
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(fc+,w,Oj_) 



Figure 2. The information encoded in Q^{s, z, ij^) is exemplified here. The incoming quark creates 
a jet of invariant mass s, inside which a hadron h with momentum fraction z = pj^ /uj is produced. 
Initially, the large parton virtualities yield emissions at wider angles. This depends on s and z, and can 
be described perturbatively by Jij(s, z, fi). At smaller parton virtualities, the emission is at smaller 
angles and essentially only affects z. Here the effect of hadronization also becomes important, and 
this is described by the standard fragmentation functions D'j{z,fi). 



The dependence on z and z' in Jij is only through their ratio, since these coefficients can 
only depend on the perturbative variables associated with the partons i and j, and not on the 
hadron h. Eq. (2.30) is analogous to the matching of beam functions onto parton distribution 
functions performed in refs. [32, 33]. 

From the tree-level results in eqs. (2.15) and (2.23), we find by using eq. (2.30) that 



J-^f is, z, tij) = 2(2^)3 5,j 5{s) 5(1 - z) 



(2.31) 



This can simply be understood that at tree-level the hadron directly fragments from the 
parton i, without emitting radiation that would build up the jet. 

The main purpose of this paper is to calculate J'ij{s, z, fij) at one- loop. This completes 
the picture detailed in ref. [8] with the information necessary to relate the factorization 
theorems for semi-inclusive processes, where the jet invariant mass is probed, to the standard 
Df{z, fi) at NLO accuracy. We find that the one-loop matching coefficients J'^p are given by 



2(27r)3 



2(2^)3 



27r 

+ 8{s 

as{fJ-j)CF 
27r 



9{z)WcJ^)5il-z) + \co(^)il + z^)Coil-z) 
1 + z2)£i(l -z) + Pgg{z) lnz + e{l -z){l-z)- —5(1 - 



e{z) 



r 1 



2 '-0(7:2 



-Cn 



+ 5(s)ln(z(l-z)) 



(2.32) 

z)_ 

5{s)e{l-z)z 
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2(27r)2 



2tt 



e{z] 



+ 5{s) Ci{l-z) 



2 
2(1 



z + z 



5{l 

2\2 



1 



Z) + -^Cq 



+ Pgg{z) InZ 



(2.33) 



6(1 



Jgl\s,z,fij) _ as{fij)TF 



2(27r)3 



27r 



e{z)\\^Co(^)+6{s)ln[z{l-z)]\Pgg{z)+26{s)d{l-z)z{l-z) 



Jgl\s, Z, fXj) = J^l\s, Z, flj) , 



where the plus distributions £„ are defined in eq. (A. 2). Furthermore J'gg = Jqq and Jqg = Jqg 
by charge conjugation invariance in QCD. The coefficients Jqq, Jqqi and jTgg', where denotes 
a quark of a different flavor than g, only start at two loops. 

The are extracted from our partonic calculations of Q^^^^ i^s, z, /i) and D'^^^^\z, ^u), by 



using eq. (2.30) expanded to one loop, 
'^dz' 



gf\s,z,fij) 



k 



E 



Id/ 



.(0) 



2{27rf6ik S{s) 6(1-^) D^^^^ {z\ fij) +J^,'> ( s, ^, /ij ) 5^, 6(1- z' 



(1) 



2(2vr)^J(.)Z)p(z,/.j) + J-\^^(s,z,/ij) 



(2.34) 



Therefore, the NLO matching coefficients can be obtained through subtractions of the one- 
loop renormalized partonic fragmenting jet functions and the one-loop fragmentation func- 
tions. The calculation of these two is the subject of sections 3, 4 and appendix B. 

2.5 Relationship between Q^{s,z,ii) and the Jet Function Ji{s,^) 

In light of the discussion in ref. [8], we re-derive the relationship between the fragmenting jet 
function and the jet function, exposing all subtleties. When we sum over all possible hadrons 
h £ Tii fragmenting from a parton i and belonging to the jet, the fragmenting jet function can 
be related to the inclusive jet function Jj(s, /x), which is completely calculable in perturbation 
theory. We use the completeness relation 

r^^^ ^ ^|X/i(z))(XMz)| = ^|X,)(X,| = t, (2.35) 



where {|Xj)} is a complete set of states in the jet-like kinematic region that we are interested 
in. The factor z under the integral is needed to provide the correct symmetry factor for 
states with identical particles. This is easily seen from the following example: consider the 
case where X consists of n hadrons identical to h, i.e. X = {h{zi) . . . h(zn)}. The sum over 
X in eq. (2.35) contains a phase-space integral over the momentum fractions zi, . . . ,Zn that 
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is subject to a momentum conserving delta function and has a symmetry factor of 1/n! 

/ dzz— Y\dzi\zi,...Zn;z){zi,...Zn;z\6{l-z-y^Zi] 
Jo n\ Jo l\ \ U J 

„l n+1 / n+1 V 

/ TT dZi \zi, . . . Zn,Zn+l){zi, . ..Zn,Zn+l\^{ 1 " ) + ^) ^"■+1 

„l n+1 ^ n+1 V 

/ rr dZj |zi, . . . Z„,Zn+l)(^l, • • -^^ni^n+ll 5 ( 1 - y^^^j ) • 
■^0 i=l V / 



(n+1) 



(n+1)! 



Since we integrate over z, the hadron h is no longer distinguishable from X and should 
therefore be grouped with the rest. In the first equality, we redefined Zn+i = z and divided 
and multiplied by n + 1 to get the correct symmetry factor in front. In the final step we 
replaced (n + l)^n+i Yl^=i which is justified because all particles are identical and all 
momentum fractions are integrated over. Using the momentum conserving delta function, 
Yl'i=i Zi = I, which leads to the result. 

Applying eq. (2.35) to the fragmentation function leads to 

V /" dxxD'^{x,fi) = 1, (2.37) 
h 

which is consistent with momentum conservation and with the definition of D^{z,fi) as the 
number density of the hadron h in the parton i [5] . Similarly, in the case of the fragmenting 
jet function we obtain 

V / dzzg^{s,z,fi) = 2{27rfJi{s,fi). (2.38) 

hen, 

Combining this with eq. (2.30), leads to 

1 /"^ /"^ dx f z \ 

•^^^^'^) = 2(2^^/0 "^"'^A 

h j 

= 2(27r)3 2^ Jq ^uuj,j{s,u,n). (2.39) 

Here we introduced the variable u = z/x to disentangle the integrations and used eq. (2.37). 
This relationship with the jet function does not constrain non-perturbative physics, but pro- 
vides a cross-check of our perturbative calculation of Jij, since the quark and gluon jet 
functions are known. 

3 Quark Matching Calculation with Gluon Mass and 6 Regulator 

In this section we present the calculation of the Wilson coefficients Jqj for matching the 
quark fragmenting jet function onto fragmentation functions at NLO. These are extracted 
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using eq. (2.34), for which we need to calculate the partonic fragmentation functions and 
fragmenting jet functions at NLO. To regulate UV divergences we employ dimensional reg- 
ularization (DR) in d = 4 — 2e dimensions and renormalize according to the MS-scheme. 
Concerning the choice of the IR regulator, we note that neither an offshellness nor a fictitious 
gluon mass takes care of the IR divergences for emissions from the collinear Wilson line Wn 
in eq. (2.4). According to ref. [24], these can be regulated as follows: 



with 5 > 0. In this section we adopt a gluon mass and 6 as IR regulators. A non- vanishing 
gluon mass violates gauge invariance when the calculation involves a triple gluon vertex, 
which is however not the case here. 

A naive calculation of the graphs includes the region when i becomes soft, which was 
excluded in eq. (2.2); zero-bin contributions have to be subtracted. In these zero-bin sub- 
tractions we need to apply the 5-regulator prescription as it would have appeared in the soft 
Wilson line: 



In this framework DR does not regulate any IR singularity, and this enables us to show in a 
clean way how IR divergences get cancelled in the matching between Qq and Dq. 

In appendix B we compute the same diagrams with a quark-offshellness regulator, where 
the IR divergences from eikonal propagators are regulated by DR. The resulting J7ij turn 
out to agree with those computed in this section, as expected since these Wilson coefficients 
should be insensitive to the choice of IR regulators. After having studied in detail the IR 
structure for the case of quark fragmentation here, we perform the gluon matching calculation 
using DR for both the UV and IR in section 4. 

In our partonic calculation we replace the hadron h in the intermediate state of eqs. (2.8) 
and (2.21) by a quark or a gluon and the remainder X by the vacuum or a gluon or a quark, 
as required at one-loop order. In this section we evaluate the graphs by integrating over the 
phase-space of the parton which replaces X. In appendix B we compute the diagrams for the 
fragmenting jet function following an alternative approach based on the optical theorem. 

The Feynman diagrams contributing to and at one loop are shown in figure 3. 
By our choice of coordinates the incoming parton has no perpendicular momentum. Both 
in D and in Q the ratio of the large components of the momentum of the incoming quark 
and the outgoing identified parton is measured. For Q, in addition, the virtuality of the 
incoming quark is specified through k'^ . The graphs in fig. 3(a) and (b) correspond to a 
quark fragmenting into a quark while emitting a gluon, fig. 3(c) corresponds to a virtual 
correction and fig. 3(d) is the contribution from the quark wave-function renormalization. 
Figures 3(b) and (c) have a mirror image which we do not draw separately, they also have 
nonvanishing zero-bin contributions. The graphs for a quark fragmenting into a gluon while 
emitting a quark are obtained by interchanging the momenta p o £ in fig. 3(a) and (b). 






•-perms 
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Figure 3. Feynman graphs contributing to the partonic fragmentation function and fragmenting 
jet function that are non-zero at one-loop in Feynman gauge are shown here. Graphs (b) and (c) 
have a mirror image and (d) corresponds to the wave function renormahzation. For the partonic 
fragmentation function D]^{z) the minus component cj of the incoming momenta is fixed and the 
perpendicular components are zero by the choice of coordinates. For the fragmenting jet function the 
plus component k'^ is also fixed, determining the invariant mass of the jet. 



We work in the Feynman gauge, without any loss of generality since the sum of the 
graphs is gauge-invariant. We use QCD Feynman rules to present the calculation here and 
have verified that the SCET Feynman rules give the same result, as expected. 



3.1 Quark Fragmentation Function at NLO 

We start with graph (a) fig. 3, 
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(3.3) 



In the final step we expand in e to extract the UV divergences. 
Moving on to fig. 3(b), 
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s,pol 
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1 2 r 

x)6il -x){-+ In [x£o(l -x)- 5(1 - x) In -1 + 0(e) , (3.4) 

Ve xm^/ L a; J 



where we include a factor of 2 for the mirror graph. At the end we take the hmit 5 — t- to 
isolate the IR divergences. There we use 



..... ^ ^— = hm ; ; = £o(l - x) - 0(1 - x) In — , 

5^ol-x + 5/uj 5^0 l-x ^ ^ ^ ^ CO 



lim 



(3.5) 



which follows from the definition of Cq in eq. (A. 2), with x = x — 6/uj. 

We need to subtract the zero-bin contribution, which comes from the region where the 
gluon becomes soft, i'^ ~ X^p~ [20]. Expanding eq. (3.4) accordingly, and including the 
appropriate 5-regulator for the soft Wilson line that appears here. 
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(3.6) 



Note that p"*" = from the on-shell condition p"^ = and pj_ = 0. After expanding in e and 
separating the IR divergences by taking — 0, followed by 5 — )■ 0, we get 

^2- 



jjq{b)0 as{p)CF 



(J, bare 



TT 



5{l-x) 



r 1 



^--ln- + 21n2- + — 
e /i M 4 



(3.7) 



This result depends on the order of the limits m? — )• 0, 5 — )• 0, but the sum of the diagrams 
does not. In eq. (3.7), is absent because it does not regulate anything for this diagram 
and can simply be set to zero. 
For fig. 3(c) we find 
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(3.8) 
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After integrating the delta functions, we perform the integral by contours. The poles are 
located at 

t = F^. '* = -i^' P'f) 

which are on opposite sides of the real axis for < < uo. We pick up the first pole and 
perform the standard ^_L-integral, yielding the second last line in eq. (3.8). In the last step 
we perform the remaining £~ integral, expand in e and take the limit 5 — >■ to isolate the IR 
divergences. 

The corresponding zero-bin contribution is obtained by expanding in the region (.^ ~ 

°£=-(^)*^^-Wi0'^<'^--'^"<-) 

S 

4^ j '9^Fd{ ^) ] ^2ttY ^l- + 5){l+ + 5-iQ){l-l+ + ll-m? + \{)) 



Att J TT ' 'Jo i- + 6 J (27r)'^-2 il-m^- 5t 

q,ha,re ' 



The opposite sign of the iO prescription in + 5 — iO comes from dividing out —p < 0. This 
time the poles are at 

= i-^ = -S + iO, (3.11) 

which are on opposite sides of the real axis for i~ > 0. Picking up the second pole bring us 
to the third line, which is equal to minus the second line of eq. (3.6). Thus the zero bins from 
the real and virtual Wilson line emission cancel each other, which is no surprise because D is 
insensitive to the scale associated to the soft radiation accompanying the final parton. This 
will no longer be true for the fragmenting jet functions. 

Using on-shell wave-function renormalization for the massless quark, 



£2 _ ^2 
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2tt 2-e m2 
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from which we derive 
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The contribution to D from wave-function renormalization is therefore given by 
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(3.14) 



We will now calculate the diagrams for Dq, where the momentum fraction of the gluon 
is measured. We can obtain these from the corresponding expressions for the quark case in 
eqs. (3.3) and (3.4) by taking x — )• 1 — x. To see how this comes about, we explicitly include 
the on-shell condition for p and the definition of x in eq. (3.3) and eq. (3.4) as follows: 



J d'^-V = j d'^p(5(/)0(p°) 6(^x 



P_ 



(3.15) 



As a consequence, the quark and the gluon in fig. 3(a) and (b) are on completely equal footing, 
except that the momentum fraction of the quark was measured. Measuring the momentum 
fraction of the gluon therefore amounts to x — )• 1 — x. 
For the first diagram we obtain 
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(3.17) 



where in the last step we dropped both the plus prescription and the 5{x) terms since x > 0. 
The zero bin only contributes at x = and therefore does not need to be taken into account. 
Adding up our results yields 



D 



9(1) 
5, bare 



D 



9(1) 
(J, bare 



D 



q(a) 



g,bare 

as{n)CF 
2tt 



_|_ / jjlib) _ jjq{b)0 
' ' g,bare ij'jbare 



q,ha.re 



e{x) 



1 



+ In 



F)9{a) I ry9(b) 
g,barc g,barc 



2tt 



e{x) 



- + In — 



)'?{c)0 
g,bare 

'TT^ 

3 



q{d) 
g,bare 



lnx)p,,(x) - - -)<5(1 - x) - 20(1 - x)(l - x) 



ln(l - x)^Pgq{x) - 261(1 - x)x 



(3.18) 
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where the sphtting functions Pig are given in eq. (2.19). Note that 1/e-poles multiplying IR 
regulators cancel in the sum of the diagrams, as must be the case. In eq. (3.18) S is absent 
since it regulates those IR divergences which cancel between the real and virtual emission 
diagrams for the fragmentation function. 

3.2 Quark Fragmenting Jet Function at NLO 

Compared to the Dj case, the partonic fragmenting jet function calculation involves an ad- 
ditional 5[uj{k~^ — — p~^)], where = p]_/p~ = —p\/p~- It would seem that this fixes all 
the momentum components in real radiation graphs 3(a) and 3(b), however, as we will see, 
this is not true for the zero-bin diagram associated with (b). Starting with figure 3(a), 
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Here we rewrote 
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(3.20) 



Since all the dependence on £± is in terms of £'j_, we used spherical coordinates to perform 
the £_L-integral 



(2vr) 



d-2 



(3.21) 



which in combination with eq. (3.20) demands that s — m?/{l — z) > for a non-vanishing 
contribution; hence the factor of 6{s — m?/ (1 — z)) in the second to last line of eq. (3.19). The 
graph is UV finite, so the expansion in e is trivial. In the final step we take the limit — 
to separate the IR divergences using eqs. (A. 2) and (A. 4). 
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The calculation of fig. 3(b) combines steps from the corresponding fragmentation function 
graph and the previous diagram, 
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We take the limit rm? — t- followed by (5 — )• 0. Since the gluon mass occurs in the combination 
m? /{\ — z), different order of limits give different results for individual graphs, but the sum 
of all graphs is independent of the order of limits. We separately studied the cases z < 1 and 
s > 0, and determined the coefficient of the remaining 6{s)5{l — z) term by integration. 
The zero bin corresponding to eq. (3.22) is given by 
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Since all divergences are regulated by 5, we can simply set m? = 0, which is consistent with 
the order of limits used previously. In taking (5 — >• we used 
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which fohows from eqs. (A. 2) and (A. 4), where s = s + 6uj. 

The virtual emission graph and the wave-function renormahzation contribution are un- 
affected by the restriction on the real radiation given by the 5-function involving 
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with (r) = (c),(c)0, (d). 
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As for the fragmentation function, expressions for the diagrams where the momentum 
fraction of the gluon instead of the quark is measured can be readily obtained by z — )• 1 — z, 
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Adding up all the diagrams, 
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All 1/e-poles here are of UV origin. The 5-regulator disappears from the final result. This 
is expected, since the fragmentation function in eq. (3.18) contains no dependence on the 
(5- regulator either, and the IR divergences have to cancel in the matching between Qi and Di. 
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3.3 Renormalization and Matching 

According to eq. (2.17) the relevant countertems are 



(3.28) 



which lead to the same anomalous dimensions as in eq. (2.18). The renormalized one-loop 
quark fragmentation functions are given by 
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Our one-loop result in eq. (3.27) agrees with the fact that the /x-dependence of Qi{s, z, ji) 
is the same as in the jet function Ji(s,fi), at any order in Og. In particular, from eq. (2.29) 
we find 
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which yields the known one-loop anomalous dimension for Jq in eq. (C.8). Furthermore, Gq 
is UV finite and therefore the renormalization does not mix quark and gluon fragmenting jet 
functions. Moreover, the UV divergences in eq. (3.27) are multiplied by 5{1 — z) implying 
that the renormalization does not affect the z-dependence. 

For the renormalized one-loop quark fragmenting jet function we obtain 
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Inserting our results from eqs. (3.29) and (3.31) into eq. (2.34), we produce the matching 
coefficients in eq. (2.32). The IR divergences cancel in the matching, as they should. 

A final cross-check on our calculation is provided by the relationship in eq. (2.39) to 
0{as). We correctly find that 
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where Jg^^(s, //) is the one-loop piece of the renormalized quark jet function in eq. (C.4). 
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Figure 4. Feynman graphs contributing to the gluon fragmentation function and the gluon frag- 
menting jet function at one-loop. We have not shown virtual diagrams here as they are scaleless and 
therefore trivially vanish in DR. Graphs (a) through (c) correspond to Gg and graph (d) to Gg''^^ ■ 
Graphs (b) has a mirror image. 



4 Gluon Matching Calculation in Pure Dimensional Regularization 

In this section we calculate the matching coefficients J^gj, using dimensional regularization 
for both UV and IR divergences. Since we cannot distinguish UV and IR divergences, we will 
not be able to check that Dg and Qg have the anomalous dimensions in eqs. (2.18) and (2.28) 
and that the IR divergences cancel in the matching eq. (2.34). We have verified both of these 
statements for the quark case in section 3, where we used dimensional regularization for the 
UV and a gluon mass plus a 5-regulator for the IR. Here, for simplicity, we will assume the 
anomalous dimensions in eqs. (2.18) and (2.28) and verify that the IR divergences cancel in 
the matching. An equivalent procedure would be to assume that the IR divergences cancel 
in the matching and verify that we obtain the expected anomalous dimension for Qg. An 
additional check on the UV-finite part of our calculation will come from the relation with the 
jet function in eq. (2.39). 

The partonic graphs for the gluon fragmentation function and gluon fragmenting jet 
function are shown in figure 4. In pure dimensional regularization all integrals for the bare 
fragmentation function are scaleless, therefore they vanish. Inserting the known one-loop 
anomalous dimensions from eq. (2.18) into eq. (2.17), we find that the renormalized fragmen- 
tation functions up to one-loop are given by 
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where the 1/e-poles are IR divergences. 

For the fragmenting jet function the real emission graphs can give a non-zero contribution, 
because the measurement of k'^ now provides the Lorentz invariant quantity s = uik^ as a 
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scale in the calculation. The virtual graphs are still scaleless because = = 0. It is 
easiest to calculate using the sum over physical polarizations. Then the real emission graphs 
only contribute to physical degrees of freedom in the final state. The physical polarization 
sum in light-cone coordinates reads 
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which gives a vanishing contribution for diagrams 4(b) and 4(c) on contracting Lorentz indices 
with the 2-gluon vertex from the operator insertion; only graph 4(a) contributes to Qg. For 
figure 4(a) we find 
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and 



p /w as always. In the first line we already used 



where s _ _ 

p±+i± = to simplify the Feynman rules for the B'^ operator. In the first step we work out 
the rather tedious contractions of Lorentz indices, using eq. (4.2) and set H.^ = {1 — z)uj owing 
to 5{oj — p~ — In the second step we perform the straightforward 5- function integrals 
that are much alike quark fragmenting jet function calculation. In the final step we expand 
in e, dropping 0{e) terms. 

Finally, we calculate the mixing graph in fig. 4(d) 
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which we notice is invariant under the transformation z — )■ 1 — z. This is expected because 
the diagram is symmetric under the exchange of quark and anti-quark hues. Therefore 
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Using the known anomalous dimension of the gluon jet function from eq. (2.28), we obtain 
the renormahzed one- loop fragmenting jet function with the aid of eq. (2.29), 
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Here the 1/e-poles are of IR origin. By subtracting the one-loop fragmentation functions 
of eq. (4.1), we arrive at the matching coefficients given in eq. (2.33). Note that the IR 
divergences again cancel, as they should. 

As a final check, we verified that these results satisfy the relationship in eq. (2.39) with 
the jet function. 
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Here, we included the factor nj to account for the number of light quark and antiquark flavors. 
The Jg^\s, /x) denotes the one-loop terms of the renormahzed gluon jet function in eq. (C.4). 



5 Numerical Analysis 

Here we present plots of quark and gluon fragmenting jet functions for vr "'"-production (sec- 
tion 5.1), and a numerical study up to NNLL accuracy of single vr "'"-fragmentation in e^e~ 
collisions where a cut on thrust is imposed (section 5.2). 

As input we use the HKNS fragmentation functions [34]. For simplicity we will always 
utilize DJ* {z) at NLO, even though our formal counting (cf. table 1 on page 31) would imply 
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Figure 5. Fragmenting jet functions for the 7r+-fragnientation from a w-quark (top row) and a d-quark 
(bottom row). Shown are the LO and NLO results with the corresponding perturbative uncertainties, 
as explained in the text. The left panels display zGi {s"^^^ , z , n) for s'^"* = (3GeV)^ as function of z, 
at the jet scale ~ v^s™*. The right panels show the same curves relative to the LO. 

the use of LO fragmentation functions at LO, LL and NLL. Consistently with ref. [34], we 
set as{mz) = 0.125, used two-loop running, and matched continuously across the b- and 
c-quark thresholds at mf, = 4.3 GeV and rric = 1.43 GeV. In our analysis we will not include 
the effects of the uncertainties associated with the fragmentation functions or as{mz)- 

5.1 Fragmenting Jet Functions up to NLO 

In this section, we will show plots for the fragmenting jet functions, using eq. (1.5) and our 
results for Jij. We will study the (dimensionless) integral of the fragmenting jet functions 
over the jet invariant mass up to ■s'^^* 

1 r""' 

gf(.-*,z,/.) = 2^y^ dsgHs,z,fi) (5.1) 

for scales n ~ where there are no large logarithms involving the variable s. Our 

plots show the effects of constraining the invariant mass of the jet both on the normalization 
and on the shape of the z-dependence. In figure 5 we show the vr "'"-fragmenting jet function 
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Figure 6. Same as figure 5 but for rf-quarlc (top row) and giuon (bottom row) initiated fragmentation 
of 7r+ witli ~ \/s™t = 10 GeV. 

Qf^{s™^, z, fi) for a representative fixed value of s'^"* = (3 GeV)^ as a function of z. This clioice 
of cut is motivated by the next section, where we consider e~^e~ colhsions at Q = 10.6 GeV. 
The top row in figure 5 corresponds to fragmentation from a n-quark, and the bottom row 
to fragmentation from a d-quark. The u- and (i-quark results are identical to those for the d- 
and n-quark, as a consequence of the fit ansatz in the extraction of the HKNS fragmentation 
functions. The left panel shows z Qf'^{s'^^^ , z, fi) and the right panel displays the corrections 
relative to the LO result ^^^\s^^^ , z , fi) = DJ^ {z, fi) 6{s™^). To obtain the fragmenting 
jet function at a scale different from the jet scale, one can use the ROE in eq. (2.25), which 
does not affect the z-dependence. 

The plots are cut off at z = 0.2, because our factorization formula is not valid for too 
small values of z. For 1 — z <C 1, eq. (5.1) contains large logarithms oi 1 — z [see eq. (2.32)], 
which should be summed. We avoid this issue by cutting off the relative plots at z = 0.8. As 
is clear from the absolute plots, z ^]^^(s™*, z, /i) is very small in this region anyways. 

To estimate the uncertainty from higher-order corrections, we vary /x between y/sjl and 
2y^, which is shown by the uncertainty bands. At LO the scale variation is simply that of 
the fragmentation function, whose maximum and minimum are obtained at ;U = -^/2,2y^ 
with central value at /x = yfs. At NLO, the maximum and minimum for the /x variation do 
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not occur at these values, due to the double logarithms in the Jij. We therefore sample over 
y/s/2 < /i < 2-y/s to determine the uncertainty band, and take the central value to be the 
average of the maximum and minimum variation. 

As can be seen in figure 5, for the u-quark the uncertainty at NLO is less than at LO 
and the uncertainty bands overlap, so perturbation theory is well-behaved. For the d-quark 
the uncertainty bands overlap as well, but they become rather large for z > 0.7, both for the 
LO and NLO result. For these values of z, -DJ^(z,/i = 1 GeV) is tiny [34] and its increase 
is due to the running up to /i ~ which leads to these large uncertainties. For larger 

scales this effect decreases: in the top row of figure 6 we show the d-quark fragmenting jet 
function for ~ -y/s™* = 10 GeV, where the uncertainties are smaller and decrease from LO 
to NLO. The gluon fragmenting jet function exhibits the same feature, which is why we again 
choose // ~ ^/s^ = 10 GeV for the plots in the bottom row of figure 6. This scale choice is 
also relevant for hadron colliders, where one would expect to see gluon-initiated jets of higher 
invariant mass. 

5.2 Pragmentation at e+e Collisions, v^ith a Cut on Thrust 

We will now show results for e+e^ — t- Xt:^ in the dijet limit, where the fragmentation variable 
z is measured. We remind the reader that in our choice of frame the jet's perpendicular 
momentum vanishes. 

Including the RGE evolution kernels in eq. (1.6), the resummed cross section is given by 

dV 



H{Q'^,flH)UHiQ'^,fJ-H,fJ-) ^ 2(2^)3 J "^^^ "^^^ 

q={u,u,d,d,s,s} 

X j OqiSa " s'^, Z, ^j) [/^(s'„, flj, fl) j ds';, Jqisb " 4- f^j) f^j(4> /^J' /^) 

X ldkQSr{QT-'^^^^-k,iJs)Us{k,fis,Pi), (5.2) 

where we sum over light quark flavors. This formula only describes the singular contribution 
to the cross section, which goes like ~ (In'^ ''")/''" for small r. The nonsingular contribution is 
suppressed by 0{t) relative to the singular one and we therefore neglect it in our numerical 
analysis of the dijet limit r <C 1. The double logarithms of r are summed by evaluating 
the hard, (fragmenting) jet, and soft function at their natural scales (iih — Q, fJ-j — Vt'Q, 
~ tQ, respectively) and then evolving them to an arbitrary common scale /i. In terms of 
the Fourier conjugate variable of r, denoted by y, the cross section takes the following form 

In^ lny{aslny)'' + {oslny)^ + asioslny)'' + . . . , (5.3) 
dy 

where the k runs over the positive integers. The terms on the right-hand side correspond to 
the LL, NLL and NNLL series. For 1 — z ^ 1 we cannot trust the convergence due to the 
large double logarithms of 1 — z, as discussed in the previous section. 
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matching 




r 

J- cusp 




LO 


0-loop 






1-loop 


NLO 


1-loop 






2-loop 


LL 


0-loop 




1-loop 


2-loop 


NLL 


0-loop 


1-loop 


2-loop 


2-loop 


NNLL 


1-loop 


2-loop 


3-loop 


3-loop 



Table 1. Order counting in fixed-order and resummed perturbation theory. 

To calculate the cross section in eq. (5.2) at a specific order, we need the input summa- 
rized in table 1, where "matching" refers to the fixed-order contribution, to the non-cusp 
anomalous dimension, Tcusp to the cusp anomalous dimension and (3 to the QCD /3-function. 
The evolution factors and the one-loop hard, jet and soft function are all known and collected 
in appendix C. Our one-loop calculation of the matching coefficients Jij is the remaining 
ingredient necessary to sum the logarithms of r to NNLL order. The convolutions of plus 
distributions are carried out using the identities from appendix B of ref. [35]. 

We will now address our choice of scale for the central value of the cross section, as well 
as the scale variations used to estimate the perturbative uncertainties. We start by observing 
that the hard function for e^e~ — >• dijet is the square of a time-like form factor and contains 
large vr^-terms for = Q from ln^(— iQ/^/^-). To improve convergence we resum these 
vr^-terms by taking fin = —iQ [36-39]. 

Following ref. [40], we observe that there are three distinct kinematic regions where the 
resummation of the logarithms of r must be handled differently: 

1) /U/^~-i(5, fij ~ y^AqcdQ, ^5 = Aqcd, 

2) ^ -iQ , ^ , fJ-s ^tQ, 

3) ifiH = tJ-j = fJ-s - Q ■ 

Here we shall focus on region 2). However, our choice of scales and the scale uncertainties 
are affected because we get close to the regimes 1) and 3). In region 1) r is small and the 
soft scale becomes of order Aqcd- The factorization theorem in eq. (5.2) remains valid, but 
non-perturbative corrections to the soft function must be taken into account, which can be 
done using the methods of refs. [35, 41]. On the other hand, for large r the resummation 
becomes irrelevant (except for the tt^ resummation) and nonsingular corrections should be 
taken into account. As was observed in ref. [40], there is an important cancellation between 
the singular and nonsingular cross section in eq. (5.2) in the limit r — ?■ 0.5, which requires 
the scales to merge in region 3). 

Our choice of scales should smoothly connect these regions, which we achieve using profile 
functions. This approach has been previously used to analyze the B — )■ Xsj spectrum [35], 
thrust in e~^e~ — >• jets [40] and Higgs production through gluon fusion with a jet veto [42]. 
Here we use the same profile functions as in eqs. (2.52) and (2.53) of ref. [42], and estimate 



- 31 - 




Figure 7. The cross section for e+e^ — > Xtt^, for Q = 10.6 GeV and z = 0.6, as function of the cut 
on thrust r'^"'. The left panel shows the resummed results at LL, NLL and NNLL order. The right 
panel shows the LO and singular NLO compared to the NNLL. In the right panel we switch off the 
TT^ resummation for the NNLL to show how it merges with the singular NLO. The bands correspond 
to the perturbative uncertainties as explained above eq. (5.4). 



the perturbative uncertainty by taking the envelope of the three independent variations of 
profile parameters in eq. (2.55) thereof (with the replacement niH — >• Q)- Our central curve 
corresponds to the following choice of profile parameters: 

fi = Q, eB = es = 0, /xo = 2GeV, n = , ra = 0.25 , T3 = 0.5 . (5.4) 

We will show here plots for the cumulant of the cross section in eq. (5.2) 

where the dijet limit is imposed by requiring r < r'^"* <^ 1. We only consider the contribution 
from the light quark flavors (uds). A strong cut on thrust almost entirely removes the f)-quark 
events [6]. The c-quark contribution to tt^ fragmentation is small and therefore neglected. 

In figure 7, we show the cross section for e~^e~ Xir^ in eq. (5.5) for the Belle cm. 
energy Q = 10.6 GeV and a representative value of z = 0.6 as a function of r'^"*. We checked 
that our numerical results are not sensitive to small variations of the scale Q. The left 
panel shows the result at LL, NLL and NNLL order in resummed perturbation theory. The 
uncertainties at LL and NLL are rather large, so only at NNLL we obtain a useful prediction. 
Furthermore the LL band is smaller than the NLL one and does not overlap the NNLL, which 
indicates that the LL is not reliable. The NLL and NNLL results are compatible within their 
uncertainties. 

In the right panel of figure 7, the LO and singular NLO cross sections are plotted together 
with the NNLL without 7r^-resummation. The singular NLO is obtained from eq. (5.2) by 
setting fiH = fJ-J = IJ-S = fJ-- The remaining nonsingular terms that are present in the full 
NLO are suppressed by 0{t) relative to the singular ones. We take /i = Q for the central 
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Figure 8. The cross section for e'^e" — > Xtt^ , for Q = 10.6 GeV witli a cut on tlirust of r < t'^"* = 0.2, 
as a function of tlie momentum fraction z. In the left panel zd(T/dz(T'^"*) is plotted at LL, NLL and 
NNLL. The right panel shows the same curves and bands as a percentage relative to the NNLL. The 
bands correspond to the perturbative uncertainties, see text above eq. (5.4). 

curve and vary /i between Q/2 and 2Q to estimate the uncertainties of the LO and NLO 
results which do not turn out to be compatible. For large r the resummation is unimportant 
(except for the vr^ resummation, which we switched off in this plot) and so the NNLL merges 
with the singular NLO. However, below r ~ 0.2 the NNLL and singular NLO start to differ 
and below r ~ 0.1 this difference is no longer captured by the uncertainty bands, implying 
that resummation is necessary. 

In figure 8 we show results for the cross section for Q = 10.6 GeV and r™* = 0.2, as 
function of z. We plot the LL, NLL and NNLL results times z and the right panel shows the 
same results relative to the NNLL (our best result). As in figure 7, the LL prediction is not 
reliable (it does not overlap with the NNLL) but the NLL is. Finally, in figure 9 we compared 
LO and singular NLO to our NNLL result. This illustrates the effect of the thrust cut on the 
dependence on the measured fragmentation variable z. As in figure 7 the difference between 
the LO and singular NLO is not captured by their uncertainties. 

From the right panel of figure 7, one would expect resummation to be only marginally 
important for t^^^ = 0.2. The difference between the singular NLO and NNLL in figure 9 
is mainly due to the vr^ resummation, which is an overall factor. If one switches off the vr^ 
resummation, there is still a difference between the singular NLO and the NNLL, but for 
most of the plotted range this is within the uncertainties. 

We also study the impact of our results on the determination of the fragmentation func- 
tion parameters. For simplicity we only consider the contribution from the u-quark here. The 
HKNS parametrization is given by [34] 

(z,^ = IGeV) = , , " T -z< (1-z)^" , (5.6) 

where determines the normalization of D'^^ , and j3^^ describe its shape, and 

B is the Euler beta function. With M^^ = 0.401 ± 0.052, al^ = -0.963 ± 0.177 and 
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Figure 9. The cross section for e+e~ — ^ Xtt^ , for Q — 10.6 GeV witlr a cut on tlirust of r < r'^"* — 0.2, 
as a function of the momentum fraction z. Here we compare the NNLL resuh with the corresponding 
LO and singular NLO. The inset shows the same curves and bands as a percentage relative to the LO. 
The bands correspond to the perturbative uncertainties. 

(3Z = 1.370 ± 0.144 we reproduce our NNLL result that we will treat as "data". We then 
fit for the three parameters to these "data" using the LO formula for the cross section, 
da^^/dz = a^DZ^{z,fi = Q). We find that al^ and Pl"" change by ~ 30% and M^'^ 
changes by ~ 70%. This clearly shows that if we use the LO result, rather than NNLL (or 
NLO), to extract the fragmentation function parameters in the presence of a cut on thrust, 
they may differ significantly from their true values. 



6 Conclusions 

In this paper we have calculated the matching coefficients Jij{s, z/z' , fi) at one-loop, which are 
an important ingredient for factorization theorems that describe spin-averaged fragmentation 
of a light hadron h fragmenting from a light quark or a gluon i within a jet with constrained 
invariant mass. These matching coefficients contain the short-distance physics relating the 
fragmenting jet functions (s, z, fi) - that depend both on the fragmentation variable z and 
on the invariant mass s of the jet - to the standard fragmentation functions D^{z' via a 
convolution in z' . We have presented our calculation for J7ij in great detail, using various IR 
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regulators for the partonic and Dj, exposing the structure of the zero-bin subtractions. 
A powerful cross-check on our results is provided by the relationship between J'ij and the 
leading jet function, which we have derived here in detail. 

We have applied our results to study fragmentation of a vr"'" in e+e" collisions, where 
we restrict to the dijet limit by a cut on thrust. Our calculation of J7ij enables us to resum 
the logarithms of r™* = 1 — T'^"* up to NNLL accuracy in the cross-section d(T/d2;(T™*). We 
analyzed this cross section for r'^"* < 0.3 and cm. energy equal to 10.6 GeV, as in the study 
of light quark fragmentation in i?-factories on the T(45) resonance [6]. Here the convergence 
of resummed perturbation theory is better than that of fixed-order perturbation theory, and 
the perturbative uncertainties become reasonably small at NNLL accuracy. The NNLL cross 
section deviates from that at NLO for r < 0.2, and below r < 0.1 resumming the logarithms of 
r is necessary for a reliable prediction. Consistent with this observation, we have also shown 
that using cross sections at LO instead of NNLL (or NLO) could have a sizeable impact on 
the extracted numerical values of the model parameters for from fits to experimental data 
for e~^e~ — )• dijet + h. 

We leave for future work the inclusion of nonsingular terms in the thrust distribution, 
and the effects of the uncertainties associated with the fragmentation functions and as{mz)- 
This would provide a reliable theoretical framework to constrain fragmentation functions from 
i?-factory data where cuts on thrust are applied. However, the results presented have a more 
general applicability, and can be used to study fragmentation taking place inside any well 
separated jet. 

Note added: While we were completing this work, a paper [43] appeared which contains a 
calculation of the quark matching coefficient Jqq to one loop. We agree with this result. 
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A Plus Distributions and Identities 

The standard plus distribution for some function g{x) can be defined as 




with 




(A.l) 
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satisfying the boundary condition dx [0{x)g{x)]^ = 0. Two special cases we need are 

'9{x) In'^x" 



Cnix) 

C^{x) 



lim 

/3-5>0 



h - p) 



n+l 



9{x) 
x^-v 



lim 



1-r, 



In our calculations we will need the identity 



9{x) 



1 



<5(x) + £o(x)-e£i(x) + 0(e2) 



and the two limits 



lim 

13^0 



X 2 



vr 



- — 6{x) 





(3^0 X^ 



(A.2) 



(A.3) 



(A.4) 



Away from x = these identities are straightforward, while the behavior at x = is obtained 
by taking the integral of both sides. General relations for the rescaling and convolutions of 
Cn{x) and C^{x) can be found in App. B of ref. [35]. 



B Quark Matching Calculation with Offshellness IR Regulator 

Here we present the one-loop calculation of Dq{x,n), Dq{x,n), Gq{s,z,fi) and Gq{s,z,fj,) 
where the IR divergences are regulated through a small quark offshellness > 0. The 
real emission graphs for the fragmenting jet function are calculated using the Lehmann- 
Symanzik-Zimmermann (LSZ) reduction formulae combined with the optical theorem. We 
use dimensional regularization with d = 4 — 2e for the UV divergences and renormalize 
according to the MS scheme. The diagrams are computed in the Feynman gauge without any 
loss of generality since their sum is gauge- invariant, and we use the SCET Feynman rules. 
The zero-bin graphs vanish in this calculation. 

B.l Quark Fragmentation Function 

For the virtual emission diagrams, the state \Xh) in eq. (2.8) gets replaced by an off-shell 
quark. In the real graphs contributing to Dq, X is an on-shell gluon, and in Dq it is an off- 
shell quark. At the end of the calculation we expand in e and neglect terms of 0{e). Starting 
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with Dq, the real emission graph in figure 3a is given by 
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with the momentum fraction x = p /u. The contribution from the diagram in fig. 3b plus 
its complex conjugate, "mirror" graph, is 
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For the virtual graphs in fig. 3c and fig. 3d and their mirror ones we find: 
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using the one-loop on-shell wave function renormalization with an offshellness IR regulator. 
As we noted in section 3.1, the zero bin does not contribute, because the fragmentation 
function is insensitive to the soft region. Adding up all these graphs, we obtain the same 
renormalization factor Z^(x,^) as in eq. (3.28). For the renormalized quark fragmentation 
function we find 
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2tt 
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The one- loop Dq is given by the real emission graphs in fig. 3a and b, with the role of 
the quark and gluon interchanged. Therefore, from the results above, 
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This leads to Z^g{x,fi) as in eq. (3.28) and 
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B.2 Quark Fragmenting Jet Function Via The Optical Theorem 

We now move on to the computation of the one-loop real emission graphs corresponding 
to Gq{s,z,fj,) and Gq{s,z,iJ,). We use the LSZ reduction together with the optical theorem, 
along the lines of the calculations in the so-called cut vertex formalism in refs. [44, 45] . This 
procedure leads to the same result as directly integrating over the parton phase space, which 
has been employed everywhere else in this paper. 

Applying the LSZ formalism to the collinear matrix elements in the definition of Gq yields 

g^q{s,z) = I dSe^'^y-/' Jd'p^^^^tTY, [^{o\K,T,6o,v^Xn{y)]\Xq{pe,Pr)) 
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where we distinguished label and residual momenta in the collinear quark states and fields. 
The superscript "amp" indicates that the lines corresponding to the collinear quark fields 
.^n™^ and ^n™^ should be amputated and replaced by the associated spinors. Rq is the residue 
of the quark two-point function. In our calculations Rq = 1 because we chose to absorb the 
finite terms of the self-energy diagram into the wave function renormalization [see eq. (B.3)]. 
Following the optical theorem, we obtain 
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Applying eq. (B.8) to the contribution to the fragmenting jet function coming from figure 3a, 
leads to 
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where s > p^ > 0. Expanding in e and evaluating the imaginary part yields 
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In the last step we made use of eqs. (A. 2) and (A. 4) as p2 _i. q to isolate the IR divergences. 
For fig. 3b plus its mirror graph we find: 
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again using eqs. (A. 2) and (A. 4) to take the p2 q limit. 

The contribution of the virtual diagrams can be derived directly from the previous cal- 
culation of the partonic fragmentation function since 
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2(27r)3 
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The zero bin vanishes for this choice of IR regulator (1/euv ~ I/^IR = 0) but still contributes 
to the fragmenting jet function turning the IR divergences into UV divergences. 
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Furthermore, as for Dq at one loop, we obtain: 
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The UV divergences only occur in ^arc already pointed out in eq. (2.24), and we find 

that the corresponding renormalization factor Zg coincides with the one derived in eq. (3.30). 
The renormalized partonic fragmenting jet functions are then given by 
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By applying eq. (2.34), we find that the matching coefficients Jql\s,z,^) and Jql\s,z,p) 
obtained from this calculation agree with the ones given in eq. (2.32). This had to be the 
case since the Jij are insensitive to the choice of IR regulators. 



C Perturbative Results 
C.l Fixed- Order Results 

The Born cross section cJq in eq. (1.6) is given by (see e.g. appendix A of ref. [40]) 



3Q2 
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where q denotes the (anti)quark flavor, Qq is the quark charge in units of |e|, Vq^e and a^^e 
are the vector and axial couplings of the (anti)quark q and the electron to the Z as e.g. in 
eq.(A3) of ref. [40]. Here mz and Tz denote the mass and the width of the Z boson. 

The hard function for thrust, at leading order in the electroweak interactions, is the 
square of the Wilson coefficient in the matching of the quark current from QCD onto SCET, 



H{Q\iiH) = \C{Q\m)\^ . 
The SCET matching was computed at one-loop in refs. [46, 47], yielding 
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The one-loop quark jet function [48] and one-loop gluon jet function [49, 50] are given by 
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The one-loop perturbative soft function for thrust can be obtained from refs. [16, 51]: 
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The one-loop Wilson coefficients J'ij{s, z, fi), for matching the fragmenting jet functions 
onto fragmentation functions in eq. (2.30), are given in eq. (2.32) for i = q and eq. (2.33) for 
i = g. 

C.2 Renormalization Group Evolution 

The RGE and anomalous dimension for the hard Wilson coefficient in eq. (C.3) are [46, 47] 
d —Q^ — iO 
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The coefficients of the a^-expansion of rcusp(as) and ^1j{as) are given below in eqs. (C.15) 
and (C.17). By solving the RGE in eq. (C.6) we obtain the evolution of the hard function: 
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where K^{fio,fi), r7p(/Uo,/u) and are given below in eq. (C.12). 
The jet function RGE and anomalous dimension are 
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where the index i = {q,g} is not summed over. Its solution is given by [35, 51-53] 
Ji{s,n)= / ds'C/j(s-s',/io,/u) Jj(s',/io), 
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According to eqs. (2.25) and (2.28), we obtain the RGE of the fragmenting jet function as 
weh as its solution by simply replacing Ji(s, /i) — )• (s, z, /u) in the expressions above. 
The RGE of the thrust soft function is given by 
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The functions Kp(/io,/u), r/p(/io,^), K^{fio,fi) in the above RGE solutions are defined as 
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Expanding the /3-function and the anomalous dimensions in powers of a^, 
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their explicit expressions at NNLL are 
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Here r = as{^i) / as{no) and we have suppressed the superscript i on K^-, r]\. and F^. Note 
that the expressions in eq. (C.14) cannot be used across quark thresholds, where nj changes. 

Up to three loops, the coefficients of the /3-function [54, 55] and cusp anomalous dimen- 
sion [31, 56] in MS are 
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(known to hold for n < 2). 
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The MS anomalous dimension for the hard function can be obtained [57, 58] from the IR 
divergences of the on-shell massless quark form factor, which is known to three loops [59]. At 
the order we are working we only need the two-loop result. 
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The anomalous dimension of the fragmenting jet function Qi and jet function Jj are equal, so 
in particular 7g;(as) = 7j(as). These anomalous dimensions were extracted for the quark jet 
function in ref. [58] from ref. [56], and for the gluon jet function in ref. [50] from ref. [60], at 
three loop order. We only need the quark and gluon jet function anomalous dimension up to 
two loops, which are given by 
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The consistency of the RGE for the factorization theorem in eq. (5.2) implies that 75- (a^) 
27g(as), fixing the non-cusp anomalous dimension of the soft function. 
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